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niK ANALYSIS OF GEOS-3 ALTIMETER DATA 
IN THE TASMAN AND CORAL SEAS 

R, S, Mather* 

Geodynamlcs Branch 
Goddard Space Flight Center 
Grcenbelt, Maryland 20771 


ABSTRACT 

A technlc|ue has been developed for pre-processing; GEOS-3 
altimetry data to establish a model of the regional sea sur- 
face. 'Fhe algorithms, as presently used, develop models 
for a 35 x 10^ km^ area with an internal precision of tl m. 
This figure is substantially influenced by the data acquisition 
period and the sea state. There are discrepancies between 
the sea surface model so obtained and CEMO based geoid 
profiles with wavelengths of approximately 2500 km and 
amplitudes of up to 5 m in this region. TTie amplitudes are 
smaller when compared with GEMlO-based geoid determi- 
nations. However, the comparison of 11 pairs of overlapping 
passes in the region indicates altimeter resolution at the 
±25 cm level if the wavelength corresponding to the Nyquist 
frequency were 30km. In most cases, the spectral analysis 
of such comparisons indicates the existence of significant 
signal strength in the discrepancies after least squares 
fitting, with wavelengths in excess of 200 km. Regional 
studies of time varying features of the sea surface in the 
data analysis area are not currently possible due to inade- 
quate tracking support and the limited time span over which 
a dense data coverage was available. 


*()n leave of absence from the University of New South Wales, Sydney 
Australia 


FT?rT:i<T>INr, PAGE BLANK NOT Fn.M*D 
CON TEN re 

Pane 

ABSTRACT ill 

1. IN'mODl’CTION 1 

1.1 Experimental Objectives — 1 973 1 

1.2 'ITie Data Requested 1 

1.3 Hie Data Currently Available 2 

1.4 Moflification of Short Term Aims of Investigation 2 

2. REGIONAL SEA SURFACE MODELS FROM SATELLITE 

ALTIMETRY 5 

2. 1 Basic Techniques 5 

2.2 Results 7 

2.3 Conclusions 12 

3. TIME VARIA'nONS IN THE SHAPE OF THE SEA 

SURFACE 14 

4. CONCLUDING REMARKS 20 

5. ACKNOWLEDGMENTS 21 

REFERENCES 21 

ILLUSTRATIONS 

Fi>;ure Page 

1 Model of (Juasi-Statlonary Sea Surface — Tasman 

and Coral Seas SEP77 9 

2 Discrepancies Between Sea Surface Model SEP77 and 
Marsh GEMIO Detailed Gravimetric Geoid in Tasman 

and Coral Seas 10 


ILLrSTRATIONS (CnntlnucKl) 


3 

4 

5 

0 

7 

8 

A-1 
A -2 
A-3 
A-4 
A-5 


T 


! 



Discrepancies netwoon Sea Surface Model SKP77 
(Technique 1 Solution) and the GKM9 Geoid MckIcI 


(Contour Interval - 1 m) 11 

Free Air Anomalies — Australia 12 

Discrepancies Between Sea Surface Model SKP77 
( Technique 2 Solution) and the GKM9 Geoid Mcxiel 
(Contour interval - 1 m) 13 


Differential Plots of Overlapping Sea Surface and NASA 
Geoid Models Illustrating the Kffect of Krrors in the 
Gravity Field Model — Passes fi7 and 97 (Tabl'^s 3 
and 4: Row 13), Date of Acquisition — No. fi7 on Day 


185, 1975, No. 97 on Day 259, 1975 15 

Plot of Discrepancies of 3 Overlapping Passes (Pass 
Pairs 7, 8 in Table 3) In the Coral and Tasman Seas 
From GKOS-3 After Tilt ?nd Bias Corrections IG 

Power Spectra, Expressed as % Stren^jth of Si^nal(s) — 

Equation 7 - for 3 Overlapping; Passes Plus Averaj^e 

Values for 14 Pass Pairs Listed in Tables 3 and 4 18 

Heights of Local Mean Sea Levels — Western Pacific 

Coast (Australia) 2G 

South-to-North Pass Groundtracks of the G EOS-3 

Spacecraft in the Coral Sea Test Area 27 

Sea Surface Toponraphy in the Orbital Plane, 

Northeastern Australia 28 

Flow Chart Work Schedule for Proposed Sea Surface 

Topof>;raphy Determination Off Australia 29 

Distribution of Surface Gravity Data in Australia and its 

Environs 30 





THE ANALYSIS OF GEOS-3 ALTIMETEH DATA 
IN THE TASMAN AND COHAL SEAS 


1. IN'mODl’CTION 

1.1 Experimental Objectives — 1973 

The proposal submitted In early 1973, sought to determine the quasl-statlonary 
sea surface topography off north-east Australia on a differential basis usin^; the 
altimeter of the GEOS-3 spacecraft in the short pulse mode. The 1970 adjust- 
ment of the Australian levelling network Indicated the existence of an apparent 
slope of 1.7 m in the sea surface, as sampled from geodetic levelling/tide gauge 
comparisons, sea level appearing to rise in relation to the level surface towards 
the equator ( Fig. A-1). ♦ While such a rise is not unexpected from oceanographic 
considerations, its magnitude is about three times greater than that computed 
from hydrostatic considerations using temperature, pressure and salinity data. 

In principle, the GEOS-3 altlmoter provides an Independent means of resolving 
this apparent anomaly. The original research plan. Illustrated in Figure A-2, 
called for the acquisition of short pulse mode altimetry in the test area between 
10”S and 25®S, extending 500 km to sea. Laser tracking support of such altimetry 
from a site near Townsville would provide the basis for defining the radial dis- 
placement of the sea surface above the selected reference surface as illustrated 
in Figure A-3. 

It was assessed that the available surface gravity and astro-geofietic data in the 
area, on combination with satellite altimetry data in the data acquisition region 
bounded by the parallels (0^'S; (50°S) and the meridians (140"E; 180®E) could pro- 
vide a basis for computing the regional geoid on a differential bases to »30cm. 
This experiment proposal was accepted by both NASA and the Australian (Govern- 
ment's) Hesearch Grants Committee in 1974. 

1.2 The Data Requested 

The experiment plan was mexiified as it was not possible to obtain funding for 
siting a transportable laser tracking system near the test area. In addition to 
extensive coverage of the test and data acquisition areas with short pulse mode 
GEOS-3 altimetry, the supplementary types of data sought for the investigation 
are shown in Figure A-4. 


*Note I-'igures A-l to A-7 not included in GFOS-3 final rc|iort 
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The sea surface profiles obtained from the altimetry were to be corrected for 
the effects of tides, temperature and salinity variations to obtain the quasi- 
stationary sea surface ( FIk. A-3). 'Hie surface gravity data, on combination 
with the astro-neodetlc and altimetry data in the region, was to be used to com- 
pute a differential geoid with a precision of tSOcm in the test area. 'ITils, in 
turn, required that all surface gravity anomaly data were free from systematic 
error to ^0. 2mGal through all wavelengths greater than that sought in the sea 
surface topography. An additional requirement was the control of differential 
errors in the global gravity field model used, over the test area to less than 
tO. 2mGal through wavelengths greater than GOOOkm. 

Such an approach assumes that the sea surface slope sought in the test area had 
linear gradients. 'ITiis assumi»tion was considered valid as the objective of the 
experiment was to establish whether or not a uniform slope existed in the sea 
surface as obtained in Figure A-1. 

1.3 The Data Currently Available 

iTie surface gravity data whose distribution is shown in Figure A-5, were made 
available by the Bureau of Mineral Resources, Geology & Geophysics, Canberra 
for this investigation. This data was supplemented by the latest Goddard Earth 
Models (GEM) from the Geodynamics Branch at Goddard Space Flight Center 
(Lerch, et al. 1977). A global file of 1® x 1* Iree air anomalies was made avail- 
able for this investigation by the Defense Mapping Agency Aerospace Center 
(DMAAC). 'ITie available astro-geodetic deflections cf the vertical were pro- 
vided by the Division of National Mapping, Canberra. 

Tracking data in support of the altimeter data acquisition was provided by 
Wallops Flight Center (WFC). Only C-Band, S-Band and Doppler data has been 
provided to date. A Smithsonian Astrophysical Observatory (SAO) laser tracking 
system commenced operation at an Orroral Valley site near Canberra on the 
south east coast of Australia in mld-197G. No tracking data in support of altim- 
eter data acquisition in the test area has been received to date. 

iTie first passes of altimeter data in the data acquisition region were received 
in August 1970. 'ilie number of passes received by February 1977 was 44. This 
number had increased to 107 by September 1977, covering the period from launch 
until April 1970. The details of the altimetry data used in this study are summar- 
ized in Table 1. 

1 . 4 Modification of Short Te r m Aims of Investigation 

It was Intended to determine the differential quasi-stationary sea surface top- 
ography using the relation 


2 


f, - h - N 



T 


T 

I 

I 


T 



( 1 ) 


where h is the height of the stationary sea surface and N that of the K<-*old above 
the selected reference surface, the former belnn determined from the altimetry 
and precise regional tracking data, while the latter is obtained from a combi- 
nation of the available surface (gravity anomalies, astro-neodetic deflections and 
the altimetry data in the data acquisition region, A careful study of the problem 
(Mather 1975; Mather, et al. 1970) called for a revision of this procedure for 
the following reasons; 

(a) Altimetry dat a were suljject to o rbit al e rrors radially. 'Hie orbits implied 
from the data on the altimetry data tape in the test area were subject to 
radial error which could be many tens of meters (e. , Mather, et al. 1977, 
p. 30). 

(b‘ Inadequ ate surface gravity coverage. The surface gravity anomaly data 
needs to be carefully controlled by a standardization network at least as 
;{ood as IGSN 71 (Morelli, et al. 1971) if it were to play a meaningful role 
in determinations of uniform gradients in the sea surface topography in the 
test area, 'ilio gravity data should also extend over the entire data analysis 
area in order that an adequate coverage were available for a diffeiential 
geold determination in the test area. 

(c) All data is re lated to t he sea surface and not the geold . Oceanographic 

evidence for the discrepancy between the sea surface and the geoid indicate 
magnitudes of up to ±1 1/2 m. Most of this discrepancy (over 107 ( of the 
power) appears to have the characteristics of a second degree zonal har- 
monic (Mather 1975, p. 67). If this were established to be the case, it may 
be possible to reduce the magnitude of the sea surface topography by solving 
for a differential model. Low degree harmonics in the quasi-statlonary sea 
surface topography can be obtained directly from altimeter orbit analysis as 
described in (Mather, et al. 1976a). 'Hie techniques proposed in this paper 
can be implemented without making any assumptions about the nature of the 
sea surface topography. 

All surface gravity anomaly data currently available are flawed in the context 
of geoid computations due to regional elevation datums not necessarily coinciding 
with the geoid with a precision better than il m. In a regional study of the type 
originally envisaged, all the land gravity anomaly data on the Australian con- 
tinent are controlled by the Australian Gravity Standardization Network (ANGN) 
and the Australian Height Datum (AIID) (Mather, et al. 1976b). 'Fhe latter is not 
a freely adjusted level network, being distorted to fit local sea level. 'Flie re- 
sulting gravity anomaly data bank is subject to long wave errors with wavelengths 
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up to 5000km and amplitudes of up to l/2mGal. A new Kravity anomaly data 
i)ank for sea surface topography studies (Al'SCAl) 70) was prepared with a view 
to mlniml/.e the effect of such errors, the rcsultinu data set i)cinK related to the 
.lervls liay datum level surface, l-'actors taken Into account in the preparation 
of Al’SCiAD 70 are the followinj^: 

• the even decree harmonic effectc caused by usinu free air anomalies In lieu 
of gravity anomalies; 

• the effect of the atmosphere; an<l 

• the non-RcocentrIcity of Reodetic coordinates used In computinR normal 
Rravity. 

'rhe data set Al'SGAI) 70 was assessed as beinR free from errors with wave- 
lenRths Rreater than 5000 km and amplitudes In excess of 0.2mGal and therefore 
adequate for studyinR linear sea surface slopes in the test area. 

The same dcRree of certainty does not extend to the quality of the Rravi*y data 
in the ocean areas fallinR with the data analysis rcRion. While documentation 
has still to be produced rcRardlnR its quality, it Is commonly held that its pre- 
cision Is at least one order of maRnitude inferior to that of land Rravity data, 
with a slronR possibility of the errors havinR siRnlflcant wavelenRth. The use 
of such data in sea surface slope determinations may well produce distortinR 
effects in excess of t50 cm on ocean Reold computations. Consequently, it was 
considered necessary to formulate the solution to this problem takinR all the 
above factor;., into consideration as outlined in the appendix. 

'Hie determination of a linear sea surface slope over a 2500 km distance can be 
solved If refined orliits (in this case, differential radial errors over the test 
area of t50cm) and a Rlobal Rravity field model of equivalent precision throuRh 
wavelenRths Rreater than 5000 km were available. Neither of these types of data 
are available in the test area at the present time. Attempts are still beinR made 
to resolve this problem. 

In the interim, the altimeter data provided to date is beinR analyzed to study the 
followinR problems: 

(i) Determination of the shape of the sea surface on a rcRional basis from the 
satellite altimeter data, usInR the information provided on the altimeter 
data tape. 
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(ii) 'Hie cohcrencp of sea surface proflloh by comparing overlapplnu passes after 
tnte>?ral multiples of I>30 revolutions of GKOS-3, and sbidyinK the power 
spectrum of the discrepancies. 

The study at (I) is reported in Section 2 and that at (it) in Section 3. 

I 

2. REGIONAL SEA SURFACE MODELS FROM SATELUTE ALTIME HtY 
2. 1 Basic Techniques 

'rhe orbits implied from the altimetry data tapes are subject to radial errors 
varvlnn from a few meters in most cases, t«) in excess of 700m in one instance 
(Mather, et al. 1977, p. 30). This latter pass overlapped another pass 37.18 
days earlier and noi more than 5 km away. While the sea surface heights were 
discrepant at the 700m level, the fit of one pass to the other with allowance for 
corrections for tilt (n) and bias (b) by least squares ^ave ros'duals which had an 
rms of H51 cm (ibid., p. 34). 

’ The inspection of other pairs of overlapping passes showed that a basis existed 

for determining a regional mo<lel of the sea surface with a resolution of at least i 

: tl m from the altimetry data tape usinu the following assumptions: ' 

(1) Orbital errors greater than tl m can be adequately modelled by corrections ■ 

b for bias and c for tilt. 

(li) The sea surface was radially stationary during the perioti of data acquisition. 

'rhe maximum pass le^nth was approximately 3500 km. 'iTie assumption at (1) 
would be questionable for such lon>j groundtracks if the effect of gravity model 
errors with shorter wavelengths were to significantly affect the radial component 
of orbital position. However, computations appear to Indicate that the contrib- 
• ution of this effect is likely to be less than t20cjn (Wagner 1977). 

i ^ 

j As the ocean tide amplitudes in deep oceans are not expected to exceed 30cm, a 

i ■ quasi-stationary differential motiel of the sea surface with a precision of tvem ! 

(l'> 30) can be obtained by adopting one of the techniques described below if the 
following assumptions were valid: 

(a) Orbital errors with wavelengths >7000 km contributed less thanf/3crrr 
to the error spectrum. 

(b) Time variations in sea surface topography were less than it/V^cm. 
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(c) V'arlations In the height over the area adopted as a Junction (cross- 

over) point were less than t;t/\/7cni. 

(d) The ocean tide amplitudes were either ellmlnatwl by modelling or, alter- 
nately, too small to affect the adjustment (e.n. , f>nly 30 cm when the pre- 
cision sought Is tSOcm). 

'Pwo techniques of adjustment su^Kcst themselves with strong analogies to con- 
ventional geodetic levelling (Mather, et al. 1977, p. 37). Roth techniques Incor- 
porate assumption (M) above. In Technique (1) , It is assumed that the Internal 
fidelity of the pass is not In question. 'Phe corrections required to fit the j-th 
pass of a network of N passes to the true sea surface arc a bias bj and a tilt 
represented by the tirade cj . If the 1-th sea surface data point recorded at time 
tjj , was fjj at the k-th junction point whose true sea surface heluht was J■^., with 
quasl-statlonary component and temporal variation Af(tjj), It follows that 
the following relation holds: 

fk ” ^ok “ fij ‘■’j (2) 

t ij belnn the time at which the first element in the pass was recorded and Atj 
the total duration of the j-th pass. In view of the uncertainties associated with 
current ocean tide models In the region, the entire network of passes was adjusted 
In the first instance using observation equations of the form 

fk = fok + "k 

Instead of the first equality at (2), Vk being treated as a normally distributed 
quantity. Thus, If the 1-th element of the j-th pass and the 1-th element of the 
m-th pass both provided estimates of f|^, a network of observation equations of 
the form 

V - bj - b,„ + CjUjj - t|j)/Atj - - t|„^ ),'At,„ + (4) 

where v Is the residual to be minimized. TTie resulting set of observation equa- 
tions are solved by least squares for the biases bj and grades Cj for all passes 
which traverse at least one junction point. By Its very nature, the resulting 
model for the quasl-statlonary sea surface Is insensitive to 

• datum; and 

• errors in orbit Integration which are factors of position but not time. 

In Technique (11) solutions, the internal fidelity of a pass Is no longer assumed. 
The basic assumption made Is that the difference in sea surface height between 
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adjacent Junction points «mi a single pass is subject only to random errors, if 
fij, were the trial sea surface heit;ht8 at two adjacent Junction points cor- 
respondiiiK to the i-th and 1-th |x)tnts «)n the respectively on the m-th |>ass at 
which the input sea surface heights were f„„, fi,,, , the resulting observation 
equation is of the form 

V - t (f,.„ (5) 

where Afj , Afj^ are the desired corrections to the trial sea surface heights 
f,j, at the j-th and k-th Junction points respectively. 

Both techniques were used to determine 

• the heljjht of the quasl-stationarv sea surface at each Junction point; and 

• the corrections for l)ias and tilt (grades) per pass needed to fit the orl)its 
the resultin;; sea surface model. 

For a discussion of the two techniques, see (ibid., p. 39). As summarized 
therein, tests sh .-.e - that the stabllitv of solutions for bias and tilt c«)rrections 
were not Jeopat jn the case of passes lonjjer than 300 km if the size of the 

Junction points (i.e., crossovers) were increased to 100km squares, in this 
case, ull sea surfac” heights within such a square were treated as estimates 
of the heltcht of a single point. The only effect was to increase the residual noise 
in the adjusted system. It is estimated that increasing; the size of a junction 
point from a 20 km square to a 100km square increases the system noise l)V 
tHO cm (ibid. , p. 40). 

2.2 Jesuits 


An initial analysis of the GFX)S-3 altimeter data in the Tasman and Coral Seas 
was performed in March 1977 (ibid.) usinn the 44 passes available at the time 
(Table 1). The internal noise in the system of observations for this solution with 
102 Junction points was tO. 8 m. The resulting quasi-statlonary sea surface 
model for the epoch April - September 1975 (MAB77) was compared with the 
Marsh-Vincent Rravimetrlc neoid based on GKMO and provided with the altimeter 
data tape. The resulting discrepancies were hl^jhly correlated with position 
(ibid. , p. 44), the sea surface model beinn systematically biased over the area 
of comparison in relation to the K<-’oid model. Similar trends are obtained when 
comparing GKM8 to GKMIO (Marsh 1977). A major discrepancy is noted in 
the region of the Lord Howe Rise in the Tasman Sea. 

A second solution was obtained for the shapt' of the sea surface from the data 
available in September 1977 and described in Table 1. The resulting sea surface 
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model - SKI’77 - shown In l''l>{ure 1, covers a area than MAH77, the 

additional ix'nion in the north east comer taking in the New Hebrides Trench 
with depths in excess of 8000 m. 

Fijture 2 illustrates the resultinti discrepancies on comparing th“ sea surface 
model MAH77 against a gravimetric computed by Marsh and based on 

GKMIO (Marsh 1977). While these discrepancies have an rms of t3. 2 m over 
the entire re>?ion without resorting to any selective elimination of data, as 
opposed to an rms of tl.Om in the case of the MAn77/Marsh-riKMfi ^('‘uid com- 
piirisons (Mather, et al. 1977, p. 4r>), the pattern of contours in both cases is 
similar in common areas. l\vo observations of sinnlficonce can be made in the 
case of the SKP77/Marsh-(iKMl 0 comparist'ns: 

(i) The GFMlO-based Kw>id has discrepancies with the sea surface model 
SKP77 with wavelengths of approximately 4000 Um. 

(ii) The largest discrepancies are correlated with sea floor features: 

- the New Hebrides trench (sea surface up to 10m lower than neold model); 
and 

- the Lord Howe rise (sea surface hi^iher tbJ'jr ;jeoid model). 

In evaluating comparisons of sea surface and ^[eold models in this region, it 
should be recognized that the regional surface gravity field is represented by 
only a sparse data set in the Tasman and Coral seas of questionable quality. The 
discrepancy patterns reflect the absence of high frequency representation of the 
surface gravity field in a r. t.ion where the gravity anomalies have a range of 
300mGal. It fcTovvs that it would be preferable to use satellite based gravity 
fields for such comparisons in poorly surveyed ocean areas, e.g. , GKM9 (Lerch, 
et al. 1977) and thereby extract sea surface features with wavelengths less than 
2000km using comparisons of the type illustrated in Figure 2. 

Figure 3 Illustrates the discrepancies between Technique 1 solution (i.e., relative 
geometry of each pass maintained with corrections only for bias and tilt) and 
GEM9. Table 2 illustrates in summary, the statistics of comparison between 
Technique 1 and Technique 2 solutions comprising the SEP77 sea surface models, 
with both the Marsh/GEMG and Marsh/GEMIO geoid models as well as GEM9. In 
general Technique 2 solutions have smaller rms residuals on comparison. This 
apparent index of quality needs closer examination. 

An examination of Figi:ro 3 indicates that, in addition to the two features apparent 
in Figure 2, significant discrepancies are also obtained along the continental shelf 
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Fi^re 2. Discrepancies Between Sea Surface Model SEP77 and Marsh GEMIO Detailed 

Gravimetric Geoid in Tasman and Coral Seas 
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margin east of New South Wales. The Kr^^vity anomaly field in this re^jlon is ab- 
normal with a north-south belt of lar^e negative fjravlty anomalies lyins immed- 
iately to the west of a similarly oriented belt of large positive anomalies along 
the margin of the narrow continental shelf (I'ig. 4). A similar study in the case 
of a Technique 2 solution ( I'ig. 5) (each pass adjusted in sections between junction 
points and not as entities) shows a tendency of the shapes of such dominant 
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l-’isure 4. Free Air Anomalies — Australia 


features to be blurred due to the excessive freedom allowed in the adjustment 
of observations. Consequently, larger features in the sea surface with wave- 
len(fths comparable with the spacing between junction points, tend to be ab- 
sorbed into the residuals as orbit errors. 


2. 3 Conclusions 


It appears that a basis exists for obtaining models of the sea surface with a 
precision equivalent to that underlying the assumption of stationarity of the sea 
surface, over extents of GOOOkm*. The correctness of such moaels in a global 
context depend on the precision with which harmonics of the gravity field with 
longer wavelengths are known. 


It is important that the internal geometry of altimeter passes be maintained 
over the test area in order that the geometry of intermediate features in the 
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sea surface are not blurred by absorption into the system errors. 'ITie precision 
achievable by from the techniques suKKe»*ted are limited by the followin^f factors: 

(i) Tidal uncertainties, not expected to exceed t30cm, 

(ii) Mesoscale variations in the sea surface topo^jraphy of up to 4:50 cm with 
decay times of 10^ days. 




(1<1) size adopter* lor junction points - a factor controlled by computer dependent 
factors, 

tiv) 'ITie extent of data coverane In the region, 

ITie minimum practicable data acquisition period In the case of GKOS-3 Is 25 
days. AssumlnK that adequate daUi coverage and tidal models are available, It 
follows that a system of 50 km junction points should be adequate for obtalnlnj? a 
t25 cm m(Klel of the sea surface for the study of time varvinj^ features. Idie un- 
certainties In the resulting model can be expected to Increase with the data 
acquisition perlotl and decrease with the extent of coverage. 


3. TIMK VAHIATIONS IN HIK SIIAPK OF THK SEA Sl’IlFACE 

WTille no results have been oljtalned to date for determinations of differential 
quasl-statlonary sea surface topojjraphy In the test area for reasons ^Iven In 
Section 1.4, It Is possible to study chanses In the shape of the sea surface alon^ 
profiles by the least squares flttlntj; of overlapping profiles. In the case of 
GEOS-3, overlaps occur every 37.13 days (l.e. , every 520 revolutions). 

Fourteen pairs of such overlapping passes occur in the 107 pass data set used to 
prepare sea surface model SEP77. ITiese pairs of passes can be fitted to each 
other by least squares usin>? corrections for bias and tilt. ITie results obtained 
are summarized In Tables 3 and 4 (Mather and Coleman 1977), The principles 
involved in comparisons of this type are Illustrated in Fl>?ure 0. 'ITie upper plot 
shows the Marsh/GEMO }rra\imetrically enhanced geoldal model together with 
two overlappint? altimetry profiles (Nos. 07 and 97), 72 days apart. Section 2 
In the lower diaf^ram Illustrates the discrepancy behveen the two sea surface 
profiles after least squares fitting. Section 1 shows the difference between sea 
surface profile 97 and the sravimetrlcally enhanced GEMO }?eoid model. ITie 
resultinK discrepancies have wavelengths of approximately 2500km and ampli- 
tudes of about 5 m. 'fhe discrepancies between the t\\’o overlapping? sea surface 
profiles have an rms residual of fc44cm (Table 3, Row 13), the discrepancies 
exceeding the spectrum of white noise only in respect of the "bin" labelled (t = 
200 km) in Table 4. 

Fi}?ure 7 illustrates the same information in the case of 3 overlapping passes in 
the test area (Nos. 57, 02, 109). Passes 62 (37.2 days later) and 109 (148.7 day 
later) were fitted by least squares to pass 57 and the residuals and plotted in 
I'ifpjre 12 with summaries in Tables 3 and 4 (Rows 7 and 8), The rms residual 
in each case is approximately tSOcm. 
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BEFORE CORRECTION FOR TILT ANO BIAS 
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the Effect of Errors in the Gravity Field Model — Passes 67 
and 97 (Tables 3 and 4: Row 131, Date of Acquisition - No. 67 on Day 
185. 1975, No. 97 on Day 259, 1975 


BIAS CORRECTION BETWEEN NASA GEOlO 
MOOEL ANO PASS 62 = -1.6m 


e B I SNOI133S UOd 31V3S 


CO 



I N0I133S HOJ 31V3S 


Plot of Discrepancies of 3 Overlapping Passes (Pass Pairs 7, 8 in Table 3) in the Coral 
and Tasman Seas from GKOS-3 After Tilt and Bias Corrections 
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The further analv’sis of profile fllscrepancTes of the type Illustrated in Fijjures 
0 and 7 for non-trlvlal si^fnal8 In the discrepancies calls for the definition of a 
model for the expected spectrum of errors (Ki^. 8). 

If the rms residual of comparisons were to, and the numi)er of frequencies used 
in the spectral analysis were N, it follows that the contribution per frequency, tc 
a white noise spectrum is ^a/\f N, assuming the spectrum of white noise to be 
flat. 


'Hie contribution K„| per "bln" in Table 4 is 



(«) 


where n is the number of frequencies included in the bln. Table 4 sets out the 
expected noise levels per bin assuming a banded white noise spectrum, together 
with the observed contributions to the discrepancies for the five dominant wave- 
length ranges. The latter are expressed as the percentage strength of signal S 
obtained according to the relation 


S = 


L 


'Mi t Hli 




X 100 


(7) 


the index i being taken over the frequencies included per bln, A| , Bi being given 
by (Mather and Coleman 1977, p. 10) 



( 8 ) 


i being the integral number of complete wavelengths in the length 1^, over which the 
comparisons are made, ds being the sampling interval, the residual v^ occurring 
at a distance s from the commencement of comparisons. 

Non-trivial amplitudes above the expected level of white noise are obtained in 
certain wavelengths (values boxed in Table 4). 'fhe significance of these results 
needs to be assessed with caution. Possible contributing factors are the 
following: 

(i) Tidal effects . Conventional tidal models in the test area (e. g. , Hendershott 
1973, p. 81) show distance between amphidromes of up to 5000 km. If this 
were the case, tidal effects should be completely absorbed in tilt and bias 
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AVERAGE FOR ALL FOURTEEN PASS PAIRS 
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Fijjure 8. Power Spectra, Kxpressed as % Strength of SignaUs) — 
Kquation 7 - for 3 Overlapping Passes Plus Average Values for 
14 Pass Pairs Listed in Tables 3 and 4 







corrections. However, these mo<leIs are not considered error-free and 
their ability to represent fine structure of the ocean tide, If any, has not 
been established. Table 4 sets out the phase difference between the two 
passes for the M 2 tide. Current ocean tide models can represent about 
two-thirds the ocean loading of Karth gravity tides (HretreKer and Mather 
1977) In this region. 

(II) Short period orbital errors Mth wavelengths less than r)000Km. Such 
errors can be expected to be of two types; 

• Krrors which are a function of position alone. 

• f'rrors which are a function of both position and time. 

The first type of error occurs when lnt/>>{ratlnK orlilts with an erroneous 
gravity field model uslnn a fixed complement of tracking stations. 'Fhe 
second type of error occurs due to a channe In the configuration of the 
tracking stations used In lntej;ratlnjf the orbits. It Is estimated that errors 
In the gravity field model affect radial orbital position throu^jh wavelengths 
which cannot be absorbed In tilt and bias corrections, wHh amplitudes of 
less than t20cm (W'anncr 1977). In a ^ross case (Tabi-. 4, How 1), where 
the bias correction was over 700m (Table 3, Row 1), slvcniflcant discrep- 
ancies occur with wavelength equal to the length of comparisons, the ampli- 
tude belns approximately 50cm. In most other cases, the bias corrections 
are less than 10 m, lendlnH: credence to the above fijpjre. It can therefore 
be concluded that radial orbital position errors with wavelengths less than 
3500 km usln^ present day gravity models are unlikely to exceed fc20cm. 

(III) Errors In the radar altimeter . At the level of precision beln^ considered, 
this Is largely an unknown quantity. Systematic errors In the altimeter 
with periods in er^css of 6-7 minutes are absorbed In tilt and bias cor- 
rections. A study of the results In Table 4 Indicate a level of white noise 
of arounu t25cm, due to a combination of the altimeter errors and the 
sea state. Noisier rms discrepancies may be attributed to distrubed sea 
state. Attempts to correlate noisy altlmetrv groundtracks in the Tasman 
Sea with ships lo^s made available by the Australian Bureau of Meteorology 
were not conclusive. WTiile there is a tendency for noisy ^roundtracks to 
occur durln;? periods if disturbed ground meteorological conditions, the 
correlations are subject to offsets. 

(Iv) Mesosca l c variations I n th e sea surface topography with time . Significant 
non-trivlal strengths of signal are obtained more frequently in the case of 
south to-north passes than for norlh-to-south passes (Table 4). The bias 
corrections in the case of rows 1 ;md 3 are large and may result In large 
discrepancies (at the ±25 cm level) due to tracking station configuration 
changes as discussed at (ii) above. 'Hiis may also be a contributory factor 


I 


In the case of comparison 0 In Table 4, No such armiment Is ol^vious In the 
case of non-trlvlal discrepancies ol)tained In the case of comparisons 2, 4 
and 5. 


I'ne use of the non-trlvlal strength of sljjnal obtained by uslnjt Kquatlons (0) to (h) 
does not provide a basis for studylnK all contributions to the spectrum of meso- 
scale variations In sea surface topography. Ocean eddies are features of Impor- 
tance In ocean dynamics and have finite structure with variations In both position 
and time. Sltpiatures of such features will appear In the distribution of residuals 
of the form 


V 




(9) 


4. CONCM’DING RKMARKS 

ilie following conclusions can be draw»n from the analysis of the GEOS-3 altimeter 

data In the Tasman and Coral Seas; 

• Short pulse mode altimetry data appears to have adequate resolution to 
delineate variable features in the sea surface with amplitudes greater than 
20cm and wavelengths Ijetvveen 10^ and 10^ km, provided the tracklnn station 
configuration remains fixed anrl rounh seas are not encountered. 

• Such resolution can be obtained without the benefit of hl^h precision tracking 
data if the orbit Integration were performed from a fixed complement of 
tracking stations, i'he adoption of such a procedure may enable fast-varying; 
features such as ocean eddies, to be tracked without the benefit of either an 
accurate global gravity field model or a jjlobal hij^h precIslC' tracking 
network. 

• i’he determination of ocean tides can only be obtained from the f^lobal con- 
sideration of GEOS-3 altimetry data if the dominant characteristics of such 
models are of lonir wavelenuth, as Implied by the occurrence of amphldromes 
about 5000 km apart In present day representations. Such a tidal analysis 
will have to I, o preceded by gravity model Improvement to tl m with deter- 
minations restricted to areas with dense enough coverage of overlapping 
altimetry (Bretreger 1970, p. 90). 

• liic non-trivial sea surface height variatlor. ? implied In Table 4 require 
further study in relation to the available ground truth. 
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• A basis exists for cietemiininK regional models of the sea surface with an 
internal precision of tlm, where 1 is the estimated departure from the 
assumption that the sea surface Is stationary over the pi*riod of data 
analysis. Such models generated in the Tasman and Coral Seas correlate 
well with both gross features in the bathymetry and the surface gravity 
anomalies. 
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Parameters Defining the Relative Fit of Fourteen Pairs of Overlapping Passes of GEOS-3 
Altimetr>' in the Tasman and Coral Seas from April to November 1975 
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igure A-1. Heights of Local Mean Sea Levels — Western 
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’igure A-2. South-to-North Pass Groundtracics of the GEOS-3 Spacecraft 
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Figure A-4. Flow Chart Work Schedule for Proposed Sea Surface Topographv Determination Off Australia 





tribiition of Surface Gra\itv Data in Australia and its Environs 
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Kljiure 1. Mcxk'I of (^uast-Stationary Sea Surface — Tasman and Coral Seas 
SKP77 

Kl>?ure 2, Discrepancies BetA\’een Sea Surface Mcxlel SKP77 and Marsh GKMIO 
Detailed Gravimetric Geoid in Tasman and Coral Seas 

Finure 3, Discrepancies Between Sea Surface Model SEP77 (Technlfjue 1 
Solution) and the GEM9 Geoid Mcxlel (Contour Interval - 1 m) 

Fi>?ure 4. Free Air Anomalies — Australia 

Flj^rc 5. Discrepancies Between Sea Surface Model SF"P77 (Technlcjue 2 
Solution) and the GEM9 Geoid Model (Contour Interval - 1 m) 

Figure 0. ’differential Plots of Overlapping Sea Surface and NASA Geoid 

Models Illustratint>; the Effect of Errors in the Gravity Field Model — 
* asses 07 and 97 (Tables 3 and 4: Row 13), Date of Acquisition — 
No. 07 on Day 185, 1975, No. 97 on Day 259, 1975 

Flfjure 7. Plot of Discrepancies of 3 Overlapping Passes (Pass Pairs 7, 8 

in Table 3) in the Coral and Tasman Seas From GEOS-3 After Tilt 
and Bias Corrections 



Fljjure 8. Power Spectra, Expressed as 9f Strength of Si(^al(s) — Equation 7 - 
for 3 Overlapping Passes Plus Average Values for 14 Pass Pairs 
Listed in Tables 3 and 4 
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Figure A-2. South-to-North Pass Groundtracks of the GEOS-3 Spacecraft in the 
Coral Sea Test Area 

Figure A-3. Sea Surface Topography in the Orbital Plane, Northeastern 
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Flow Chart Work Schedule for Proposed Sea Surface Topography 
Determination Off Australia 
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Figure A-5. Distribution of Surface Gravity Data in Australia and its Environs 
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Klj^ure A-»J. Model of the (A«aHl -Stationary Sea Surface for the Epoch (Aprll- 
Septembor 1975) In the Coral and Tasman Seas - MAR77 

Figure A-7. Discrepancies Bet^^’ecn Sea Surface Model MAR77 c-.nd the Marsh 
Gravlmetrlcally Enhanced GEMG Geold Model In t>'< Tasman and 
Coral Seas 
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